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Optical fields can induce forces between microscopic objects, thus giving rise to new structures of 
matter. We study theoretically these optical forces between two spheres, either isolated in water, or 
in presence of a flat dielectric surface. We observe different behavior in the binding force between 
particles at large and at small distances (in comparison with the wavelength) from each other. This 
is due to the great contribution of evanescent waves at short distances. We analyze how the optical 
binding depends of the size of the particles, the material composing them, the wavelength and, above 
all, on the polarization of the incident beam. We also show that depending on the polarization, the 
force between small particles at small distances changes its sign. Finally, the presence of a substrate 
surface is analyzed showing that it only slightly changes the magnitudes of the forces, but not their 
qualitative nature, except when one employs total internal reflection, case in which the particles are 
induced to move together along the surface. 
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I. INTRODUCTION 

Some time ago, it was demonstrated that optical fields 
can produce forces on neutral particles, 1,2 since then this 
mechanical action has been used in optical tweezers 3 and 
more recently in optical force microscopy, 4,5 as well as 
in manipulating molecules 6 and dielectric spheres. 7-9 In 
addition, the possibility of binding objects 10 through op- 
tical forces and thus create microstructurcs, cither in or 
off resonant conditions, 11-13 were pointed out. 

In this paper we wish to undertake a detailed study 
on optical forces on neutral particles, based on a rig- 
orous analysis that we have carried out 14-16 in a full 
three-dimensional configuration by using the coupled 
dipole method (CDM) of Purcell and Pennypacker 17 via 
the Maxwell's stress tensor. 18 Specifically, we study the 
forces induced by light between two spheres, either iso- 
lated in solution, or in the presence of a flat dielectric 
surface. We shall monitor the nature, either attractive or 
repulsive, of the light induced force between the spheres, 
according to the wavelength, polarization of the incident 
wave, and the size and composition of the spheres. 

In Section II we outline the calculation method em- 
ployed to determine the optical binding forces, then, in 
Section III we present results for spheres either isolated in 
water (III A) or suspended in this liquid in the presence 
of a flat dielectric interface (IIIB). 



II. METHOD USED FOR COMPUTING THE 
OPTICAL BINDING 

In a previous article 14 we showed the possibility to 
compute the optical forces on a sphere with the cou- 
pled dipole method (CDM). 17 For the computation of the 
optical binding between particles, we now use the same 
procedure, thus we shall next only outline the main equa- 
tions and the changes introduced in them to address the 
presence of multiple objects. 



Let K objects be above a flat dielectric surface. Each 
object is discretized into Nk subunits, with k = 1, . . . , K. 
Following the procedure of Ref. [17], the field at the 
(i, k)th subunit, namely, the ith subunit of the kth ob- 
ject, can be written as: 

K N l 

E(r*. u) = E (rf, w ) + £ £[S(r*, t\,u) 

1 = 1 j=l 

+ T(r 2 fc ,r^)]^.E(r^). (1) 

where a\ is the polarizability of the subunit, T is 
the linear response to a dipole in free space, 19 and S rep- 
resents the linear response of a dipole in the presence 
of a surface. 20,21 The value of the electric field at each 
subunit position is obtained by solving the linear system 
Eq. (1) written for all subunits, so that the size of the 
system to solve is Yl k =i ^k- Once the electric field is ob- 
tained, the component of the total averaged force on the 
(i, k)th subunit can be deduced from both the field and 
its derivative at its position r^: 15 



F tl (r 4 fc ) = (l/2)5Re^ 

v=l 

(u=l, 2, 3). 
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where u, v, stand for either x,y,z, and p(r^,w) is the 
electric dipole of the (i, k)th subunit due to the incident 
field and all the other subunits. Notice, that the deriva- 
tive of the field can be obtained from the derivative of 
Eq. (I). 14 Then the following relation can be written: 



(3) 



where F fc is the total force on the kth object due to both 
the incident field and the multiple interaction with the 
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surface and the other (K — l)th objects. If the object 
kth is a sphere small compared to the wavelength, the 
dipole approximation can be done, hence Nk = 1. We 
also remark that, in what follows, when we represent the 
normalized force, this means F/(47re |Ei| 2 ) where £o is 
the permittivity of vacuum and |Ej| 2 denotes the inten- 
sity of the incident beam. 



III. RESULTS AND DISCUSSION 



In Fig. 1 we represent the more complex geometry that 
we shall consider in this work. Two spheres (either di- 
electric or metallic) are embedded in water (e w — 1.69). 
Illumination with an incident plane wave takes place in 
the XZ plane at an angle 9 of incidence. When a dielec- 
tric flat surface at z — is used, we consider it separating 
glass (e s = 2.25) at z < from water [z > 0). 
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FIG. 1. The most complex geometry considered in this pa- 
per: two spheres of radius r on a dielectric flat surface. The 
spheres are embedded in water with e w = 1.69, and the rel- 
ative permittivity of the surface is e s = 2.25. The incident 
wave vector ko is in the XZ plane, and 9 is the angle of 
incidence. 



A. Particles in water 



In this section we do not address yet the presence of 
the surface (e s =£10 = 1.69), i.e. S(rf,r^-, w) = in 
Eq. (1), and the angle of incidence is 9 = 0°. Even in the 
absence of surface, we make reference to the polarization 
and thus we shall always use the terms p-polarization and 
s-polarization when the electric field vector is in the XZ 
plane and along the Y-axis, respectively. 
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FIG. 2. Normalized force in the X-direction on sphere B 
versus distance x between the centers of the spheres. Both 
spheres are of glass (ei = £2 = 2.25), with r = lOnm. 
The angle of incidence of the illuminating plane wave is 
9 — 0° and the wavelength A = 632. 8nm in vacuum. The 
full line corresponds to p-polarization, and the dashed line 
represents s-polarization. a) Force for short distances be- 
tween the spheres, the symbol + ( x ) corresponds to the val- 
ues from the non-retarded approximation for p-polarization 
(s-polarization). b) Force in far field, the symbol + (x) rep- 
resents the values from the non-retarded approximation in far 
field for p-polarization (s-polarization). 

We begin with the most simple case, i.e., the radius r 
of the two particles is small compared to the wavelength 
employed. As previously said, we then use the dipole 
approximation. We study, first, the case of two identi- 
cal spheres with E\ = £2 = 2.25, radius r =10nm at a 
wavelength A =632. 8nm in vacuum. Fig. 2 represents 
the force along the X-direction on the sphere B at dif- 
ferent positions of this sphere on the A-axis. The sphere 
A remaining fixed. We have only plotted the force ex- 
erted on sphere B, since by symmetry, the force along 
the A-axis on sphere A is the opposite to that on B. 
We observe two facts: first, the oscillation of the force 
when the spheres are far from each other, and, second, 
the strong force, either attractive or repulsive, when the 
spheres are very close from each other, depending on the 
polarization. For a better understanding of the physical 
process, using Eq. (1) and its derivative for two dipolar 
objects, we can analytically determine, through Eq. (2), 
the force on the spheres. Then, on using the fact that 
there is a plane wave in the Z-dircction (the incident 
wave is E 0i e lk ° z , where i = x or y depending of the po- 
larization of the incident field) the force on the second 
sphere can be written as: 



F x (r 2 ) = l/23?e (a 2 Ei(r 2 , u)a\E*{v u ^)^T* 



(4) 



where i = x for p-polarization and i — y for s- 
polarization of the incident field. Notice that to obtain 



2 



the force on sphere A, the indices 1 and 2 must be per- 
muted. But, even in this simple case, the exact analyti- 
cal solution of Eq. (4) is not easy to interpret. Hence, 
we make in Eq. (1) the approximation that the term 
T(r* : , r' ;, <jj)a l j is smaller than 1 (we will discuss further 
this approximation) . Now, if we used the hypothesis that 
the two spheres are identical (ai = a 2 ), Eq. (4) becomes : 

F x (v 2 ) = l/2K£ 0l | 2 3?e (|^*( r 2, n. u)j . (5) 

At short distance we can make the non retarded ap- 
proximation (fco = 0) and, as shown in the appendix, 
we have that F x (r 2 ) = — 3\aiEo x | 2 /° m p-polarization 
and F x (r 2 ) — (3/2)|aii?o H | 2 /a 4 in s-polarization. The 
points (with the symbols + and x) obtained with this 
approximation are shown in Fig. 2a and fit correctly 
the curves obtained without any approximation, as seen 
in this figure. Thus, they validate the approximation 
T(r^, r l j, oj)a.j <C f previously done. Only when the 
spheres are very close to each other this approximation 
slightly departs from the exact calculation due to the 
increase of the free space susceptibility. In fact, this ap- 
proximation assumes the dipole associated to the spheres 
as only due to the incident field, which is a good assump- 
tion when the polarizabilities are small, like for glass 
spheres. It is now easy to physically understand from 
Eq. (5) the reason of this either attractive or repulsive 
force. As the spheres are small, the scattering force is 
negligible 22 and thus only the gradient force remains due 
to the interaction between the dipole associated to sphere 
B, and to the variation of the field created by sphere A at 
the position of sphere B. In p-polarization, the field due 
to sphere A at the position of sphere B and the dipole of 
sphere B are in phase, hence sphere B is pushed to the 
higher intensity region, namely towards sphere A. In s- 
polarization, as the field due to sphere A at the position 
of sphere B and the dipole of sphere B are in opposition 
of phase, sphere B is pushed to the lower intensity region, 
namely, far from sphere A. One can observe a similar ef- 
fect in an atom mirror 23 , or on a small silver particle in 
an evanescent field. 16 

On the other hand, in the far field we obtain, from 
the appendix, the force upon sphere B as: F x (r 2 ) — 
lai-Ebx | 2 /Cq cos(fc a)/a 2 in p-polarization and F x (r 2 ) = 
— \a1E0 \ k$ sin(fcoa)/(2a) in s-polarization, with fco = 
The same explanation as before can be used 
for the sign of the force: following the phase relation- 
ship between the dipole and the field due to sphere A, 
the force is either positive or negative, hence, the oscilla- 
tions of the force F x take place with period X/y/e^. The 
phase difference X/(Ayfe^) which appears in the far field 
between the oscillations of s and p polarization, comes 
from the difference between the derivative in the com- 
ponents xx and yy of the free space susceptibility. We 
observe that the force in p-polarization decreases faster 
than in s-polarization, this is due to the absence of a 
propagating field along the X-axis in the far field. The 



magnitude of the force differs by a factor I0 4 between far 
field and near field. This is due to the strong interaction 
between the spheres through the evanescent waves. 

We can make an analogy in the near field with molec- 
ular physics. If we look at the dipole moment of the 
two spheres, we compare our system of forces with the 
interaction between two molecules. In p-polarization, as 
the dipole moments are aligned and antisymmetric, they 
produce an attractive force analogous to that between 
two orbitals p z , giving rise to a bonding state <r u . In 
s-polarization, the dipole moments are parallel and sym- 
metric, so we have antibonding states it*, where * means 
that the two spheres cannot be bound. 
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FIG. 3. Normalized force in the ^-direction on sphere B 
versus distance x between the centers of the spheres. Both 
spheres are of glass with r — lOnm, 6 = 0° and A = 632. 8nm 
in vacuum. The full line corresponds to p-polarization, 
whereas the dashed line represents s-polarization. 



We represent in Fig. 3 the force along the Z-direction. 
In this case the scattering force is predominant. The 
interaction between the spheres is now directly respon- 
sible for the oscillation of the force. Notice that when 
the spheres are far from each other, as the interaction 
between the spheres becomes weak when the distance in- 
creases, the force tends toward the scattering force upon 
one sphere due to the incident field. As this force is not 
responsible of optical binding, we are not going to discuss 
it further. 
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which is larger than that of the incident field, hence the 
dipoles associated to the two spheres are in phase and 
the force is attractive. Notice that the change of sign 
occurs both at the plasmon resonance (3?e(ei) w —2e w ) 
and when JJe(ei) w e w . When E\ is between these two 
values, the real part of the polarizability is negative. For 
a more complete discussion on this, one can see Ref. [16]. 
A similar reasoning can be done for s-polarization. If we 
now make the analogy previously done with the molec- 
ular orbitals, then depending on the wavelength in s- 
polarization, we shall obtain cither antibonding states 
7r* , or bonding states tt u . 

In the far field, for p-polarization, on using Eq. (4) and 
the approximation T(r^, t 1 -, w) a j <C 1, wc can write the 
force on the sphere B as 



FIG. 4. Non-retarded approximation normalized force in 
the X-direction on sphere B versus distance x between the 
centers of the spheres. The sphere A is of silver and the 
sphere B is of glass with r =10nm. A = 365nm (full line), 
388. 5nm (dotted line), and 600nm (dashed line). Without 
symbol: p-polarization; with symbol +: s-polarization. 

More interesting is the case of two different small 
spheres, one (B) being dielectric and the other (A) be- 
ing metallic (silver). The first fact easily observed from 
Eq. (4) pertains to the non retarded the deriva- 

tive of the free space susceptibility is real, the forces on 
spheres A and B are equal but of opposite sign to each 
other (this is no longer the case in the far field) . In Fig. 4 
we represent the force on sphere B at short distances 
from sphere A (in comparison with the wavelength), for 
A = 365nm, 388. 5nm, and 600nm. For p-polarization, wc 
observe that the force at the wavelength A = 365nm has 
a rather strange behavior as it is positive, and only when 
the spheres are almost in contact, this force changes and 
becomes similar to those at the other wavelengths. With 
the approximation T(r^, r^, oS)ctj <C 1, the force in the 
non retarded approximation can be written 

F x (r 2 ) = l/2 dTu{r2 d ^\ E^ a ^e( ai ). (6) 

We observe that the sign of the force depends of jfte(ai). 
When 5fte(ai) > (A = 600nm), which is the com- 
mon case, the dipole associated to the silver sphere is 
in phase with the applied field, so everything happens as 
for the dielectric sphere. Conversely, when 3?e(ai) < 
(A = 365nm) the dipole is in opposition of phase with 
the applied field, and hence the force becomes positive. 
But when the spheres are almost in contact, the approx- 
imation T(rf,r l j 7 u;)aj < 1 is no longer valid as shown 
when 3ie(o!i) = (A = 388. 5nm), case in which the force 
is not null but negative. This is due to the fact that 
the polarizability of the silver sphere is large and hence 
the approximation is no longer valid for short distances. 
Physically, this represents the contribution of the metallic 
sphere to the electric field acting on the dielectric sphere, 



F x (r 2 ) = ($te(ala 2 ) cos(fc a) 

- 3m(aqa2) sin(fc a))^ \E 0x \ 2 /a 2 (7) 

and the force on sphere A as: 

Fa;(ri) = (-3?e(a*a2) cos(fc a) 

- 3m( a ;a 2 ) sm(k a))kl\E 0x \ 2 /a 2 (8) 

As the spheres are small, we can take only the gradi- 
ent force as this is now the predominant one, then a 2 
is real. 15 Therefore, the forces on spheres A and B for 
the wavelength A = 600nm, where Qm(ai) is weak, are 
opposite to each other like for two identical spheres. But 
at A = 388. 5nm, where 3?e(ai) = 0, the forces on the two 
spheres are completely identical. 
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FIG. 5. Two glass spheres of radius r=100nm, 9 — 0°, and 
A = 632. 8nm in vacuum. The laser intensity of the incident 
light is 0.2W/^m 2 . Full line curves are for p-polarization. 
Dashed line curves correspond to s-polarization. a) Potential 
of sphere B normalized to ktT versus distance between the 
centers of the spheres. The height of the bars correspond to a 
normalized potential equal to 3. b) Force in the X-direction 
versus the distance between the spheres. 
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It should be remarked that if the laser intensity of the 
incident light is assumed to be 0.2W//zm 2 , the opti- 
cal forces for these small spheres are not strong enough 
to create an optical binding, since then the Brownian 
motion stays the dominant force. In this respect, the 
interest of the case of these small spheres is mainly the 
interpretative value it yields of the underlying physics. 
However, for larger radius in comparison to the wave- 
length, the forces become larger and so is the trapping 
potential. In Figs. 5b and 6b we plot the force along 
the X-axis for two dielectric spheres (glass) with radius 
r = lOOnm, and 200nm, respectively. We observe that 
with the intensity used previously (0.2W//im 2 ) the mag- 
nitude of the force is now enough to optically binding 
both spheres. We compute the potential energy of the 
optical trap by integration of the force (we take the po- 
tential energy null when the second sphere is at infinity) . 
As the two spheres are identical, the potential energy is 
the same for both. The efficiency of the trapping force 
requires it to be larger than the force due to the Brown- 
ian motion, hence the depth of the potential wells of the 
trap should be larger than kbT, T being the temperature 
of water, and kb the Boltzmann constant. Considering 
T =290K, then k b T = 4 x 1CT 21 J. We plot in Figs. 5a 
and 6a the potential normalized to the value kbT. We 
adopt the criterion that the trap is efficient when the po- 
tential well is larger that 3kbT. Hence the bars plotted 
at the bottom of the wells in Figs. 5a and 6a correspond 
to the value 3. We see from Fig. 5 that for p-polarization 
the trap is not feasible except when the spheres are in 
contact. For s-polarization, we have three equilibrium 
positions spaced out by one wavelength. This behavior 
is explained by the previous results on small spheres, and 
in agreement with experiments. 10 When the size of the 
sphere is close to one wavelength, we see from Fig. 6 that 
the depth of the potential well is larger than in the pre- 
vious case. In p-polarization there is no possibility to 
stick the spheres together, but now we have one poten- 
tial minimum of stable position. Thus, we observe that 
it is easier to trap particles when their radius are large, 
in agreement with the experiments of Burns et al. 10 . 
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FIG. 6. Same as Fig. 5 but with a radius r =200nm. 

Notice that the gravity force is 6.16 x 10~ 5 pN and 
4.93 x 1CP 4 pN for the spheres of radii r=100nm and 
200nm, respectively. Another force that exists between 
the spheres is the Casimir-Poldcr force, F c . To our knowl- 
edge, F c has often been studied either between two plates 
or between a sphere and a plate, 25 but it has never 
been established between two spheres. In the non re- 
tarded approximation and the dipole approximation for 
the spheres, F c is reduced to the dispersion force (Lon- 
don's force) which is inversely proportional to the sev- 
enth power of the distance between the dipoles (see for 
example Ref. [26]). Hence, only when the spheres are 
in contact, or at distances smaller than the wavelength 
from each other, this force might be of the same magni- 
tude of the optical forces. However, the fast decay of this 
force at distances larger than the wavelength prevents it 
to perturb the optical trap. 

In table 1 we give some examples on the limiting radius 
to get optical trapping for two identical spheres embed- 
ded in water (notice that for p polarization we do not 
mean optical trapping when the spheres are sticked in 
contact), i.e., the minimum radius to obtain one poten- 
tial minimum of stable position for the two spheres using 
the same criterion as before (namely, U > 3/ctT). 



spheres 


glass (A = 632. 8nm) 


silver (A = 394nm) 


silver (A = 314nm) 


limiting radius (p-pol.) 


123 nm 


33 nm 


180 nm 


limiting radius (s-pol.) 


85 nm 


21 nm 


50 nm 



TABLE I. Minimum radius to get one minimum position of the potential for two identical spheres. The following cases are 
addressed: glass sphere, silver sphere (both off and in plasmon resonance) The criterion of stability used is: the potential well 
depth must be larger that 3fcftT. 



We should remark that this is the limiting radius only As mentioned before, the table shows that the optical 
to obtain the first stable position, if we want to get more trapping is easier for s polarization. For the silver sphere, 
stable positions, as in Ref. [10], the radius must be larger. the value A = 394nm corresponds to the plasmon reso- 
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nance, and that of A = 314nm is for a wavelength out of 
resonance. At the plasmon resonance the polarizability 
is largest, so it is easier to perform optical binding at this 
wavelength. 
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FIG. 7. The sphere A is of glass and the sphere B of sil- 
ver with r =100nm, A — 388. 5nm. The laser intensity of the 
incident light is 0.2W/^im 2 . Plot of the potential normalized 
to kbT for the two spheres versus the distance between them. 
The full line is for p-polarization and dashed line corresponds 
to s-polarization. The potential of sphere A is without sym- 
bol, and the potential of sphere B is with symbol +. 

As a last instance, we now consider two spheres in 
free space, with radii r =100nm, one being of glass 
and the other of silver, illuminated by a plane wave at 
A =388. 5nm. We plot in Fig. 7 the potential energy 
for the two polarizations and the two spheres, since now 
this magnitude depends on the sphere material. We then 
observe that it is not possible to obtain a stable equi- 
librium since the potentials of the two spheres are now 
different. This result is explainable from the previous 
calculation on small spheres (of silver and glass), since 
then the forces were always opposed to each other at this 
wavelength. In fact, as the spheres are now large, the 
forces at this wavelength are not exactly in opposition, 
due to the larger scattering and absorbing force. Hence it 
is possible to obtain points where the potential of the two 
spheres is minimum. This happens when the forces on 
each sphere are the same and positive. In that case, the 
two spheres move in the direction of the positive A-axis 
while keeping constant the distance between them. 

B. Particles in water on a dielectric flat surface 

In this section we consider a flat dielectric surface upon 
which the spheres are suspended in water, as shown by 
Fig. 1. We compute the force along the A-axis on sphere 
B when both spheres are dielectric (glass), with 9 = 0° 
(Fig. 8). We now observe that both for large and small 
spheres the force has a behavior similar to that acting on 



dielectric spheres isolated in water. When the spheres are 
in contact, the force on sphere B is the same as in absence 
of interface, whereas when the spheres are far from each 
other, this force is slightly smaller than without interface. 
This means that the optical binding is more difficult to 
perform when the spheres are on a surface than when 
they are far from interfaces. Also, there is a change in 
the period of oscillation due to interaction between the 
spheres via the light reflected by the surface. However, in 
the case when one of the spheres is metallic (silver), we 
observe the same behavior of the forces when the surface 
is present as without it. Then, as previously observed 
for dielectric spheres, there only appears a shift in the 
oscillation and magnitude of the forces. 




FIG. 8. Force in the X-direction upon sphere B when 
the spheres are placed on a flat dielectric surface. 9 — 0°, 
A = 632. 8nm. The full line represents the force in the pres- 
ence of the surface, and the dashed line corresponds to the 
force computed without the interface. The curves with sym- 
bols + denote s-polarization, and those without symbol cor- 
respond to p-polarization. a) spheres of glass with r =10nm. 
b) spheres of glass with r =100nm. 

In Fig. 9 we investigate the potential and optical force 
on two glass spheres in front of a dielectric surface, illu- 
minated by total internal reflection (6 = 50°). Figure 9b 
shows that the force component along the A-axis always 
pushes the spheres in the direction of the wave vector 
component parallel to the surface. Hence, it is not pos- 
sible to obtain a stable equilibrium with the two spheres 
remaining fixed. But if we compute the potential of the 
two spheres together (Fig. 9a), we observe some minima 
indicating that the system can acquire internal equilib- 
rium, namely, the relative positions of the spheres can be 
kept fixed. Hence, when both spheres move impelled by 
the evanescent wave propagating along the surface, their 
velocity remains parallel to this surface, while the dis- 
tance between them keeps some particular values given 
by the position of the potential minima (cf. Fig. 9a). 
Notice that the force on the second sphere (in both po- 
larizations) has no oscillation, a very similar behavior 
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was observed by Okamoto et .al. 24 . The computational 
prediction of similar collective movements in systems of 
more than two spheres will involve long computing times 
of their relative positions by potential energy minimiza- 
tion. 



IV. CONCLUSION 



FIG. 9. Fig. 9 Two glass spheres of radius r=100nm in 
vacuum, in front of a flat dielectric surface. 6 = 50°, 
A = 632. 8nm. The laser intensity of the incident light is 
0.2W//im 2 . Full line curves are for p-polarization. Dashed 
line curves correspond to s-polarization. a) Potential of in- 
teraction between the two spheres, normalized to kbT, versus 
distance between the centers of the spheres, b) Force in the 
X-direction against distance between the spheres. The force 
curve on sphere A is without symbol, and the force on sphere 
B is with symbol +. 
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mutually coherent, the resulting standing wave pattern 
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APPENDIX A: DERIVATIVE OF THE FREE 
SPACE SUSCEPTIBILITY 

The derivative of the free space susceptibility used in 
this paper is: 

d m , . 6a / a -i \ 

—T xx (x,x ) = --^, (Al) 

d d 3a 

—T yy {x,x ) = —T zz {x,x )=jr, (A2) 

in the non retarded case, and 

JL Txx (x,x ,u;) = ^e ik ° a , (A3) 
ox cr 

^T yy (x,x ,w) = ^T zz (x,x ,co) = l -^-e ik ° a , (A4) 

in the far field, where a = (x — x ) and a = |a|. a; is the 
abscissa of the observation point, and x that of dipolc 
position. 
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